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Dear Sir: 



1 . I, L. Alison Mclnnes, declare and say I am a co-inventor of the claims of the above-identified 
patent application. I directed others and personally performed the research leading to the invention disclosed and 
claimed therein. 



2. I have read the Office Action dated April 24, 200 1 in this application and understand that the 
Examiner has rejected pending claims 1-12 and 25-27 on the basis that the specification is not enabling for the 
full scope of the claims. 

3 . The data presented below show that, using techniques described in the specification, at least five 
new polymorphisms, including single nucleotide polymorphisms (SNP), were identified in the narrow interval on 
chromosome 1 8p described in the application, which polymorphisms are associated with bipolar mood disorder 
(BP) , Thus, in addition to the polymorphisms already identified in the patent application, and using the guidance 
provided in the appUcation, several additional polymorphisms were identified that are associated with BP. 
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Association of Polymorphisms with BP in a Narrow Interval on Chromosome 18p 
AS Identified in the Instant Application and Corroborated by Subsequent Work 

4. The instant application provided data showing a positive LOD score for a D 1 8S59 allele with 
BP in a pedigree analysis; and gave evidence of an association of D18S59 with BP in a population study. The 
instant application further showed a positive LOD score for a D 1 8 S476 allele with BP in a pedigree analysis and 
gave evidence of association in population studies. In a subsequent study of linkage disequilibrium (LD) on 
chromosome 1 8 in a population sample of 69 BP-I patients from the Central Valley of Costa Rica (CVCR), the 
sameD18S59 allele was associated with BP-I. Escamilla et al. (1999) Am. J. Hum. Genet. 64:1670-1678; a copy 
of which is provided herewith as Exhibit 2. 

5. Further genotyping of the 69 affected individuals using four publicly available microsatellite 
markers delineated a segment of maximal LD with BP-I, covering about 33 1 Kb. Evaluation of a larger sample 
(227 patients and relatives, and 26 independent control trios) using these markers showed continuing evidence of 
LD and haplotype sharing in this sample for this region. Escamilla et al. (200 1) Am. J. Med. Genet. 105:207- 
2 1 3 ; a copy of which provided herewith as Exhibit 3 . 

6. Thus, the instant application provides evidence of association of at least two polymorphisms 
associated with BP. This association was corroborated by work published after the filing date of the instant 
application. These markers are in a narrow interval between SAVA5 and ga203 on chromosome 18p. Within 
this region, a segment of about 33 1 kb, and having maximal LD with BP, was further delineated. 

At Least Five Additional Polymorphisms Associated with BP were Found in the 
Previously Identified Narrow Interval 

7. Using techniques described in the instant apphcation, at least five additional polymorphisms 
were identified that are located within the narrow interval between SAVA5 and ga203, and that are associated 
with BP. 

8. As described in detail below, four new microsatellite markers, and 26 new single nucleotide 
polymorphisms (SNPs) were identified in the narrow interval on chromosome 1 8p. The results of LD analysis of 
these 30 new markers, as well as four previously identified microsatellite markers, are displayed in Table 1. Of 
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the 34 markers presented in Table 1, 16 showed association (X > 0) with BP in at least one of the two samples. 
The p-value for fxve of these 16 markers was < 0.01. All five of these markers (PH84, PH205, PH202, PH208, 
and TS30)hadestimatesof5L near 1.0,indicating that virtually all affected individuals had at least one copy ofthe 

associated allele. 



Methods 
Sample collection 

9. Two samples were analyzed. In one sample, the patient sample was composed of 227 CVCR 
BP-I individuals (including the set of 69 patients from EscamiUa (2001) that gave the original association 
evidence inl8p)and their availablefirstdegreerelatives(totalN=563). All affected individuals had at least two 
psycWatric hospitalizations with the firsthospitalization by ageSO.Asecond sample was comprised oftheseS^ 
individuals and a set of controls (52 unrelated parents of students recruitedfromtheUniversity of Costa Rica who 
were selected for CVCR ancestry [at least 5 out of 8 great-grandparents from the CVCR]). 

Radiation hybrid and STS-content mapping of markers within the candidate interval 

10. Genetic and physical mapping information was initially obtained from Whitehead histitute 
for Biomedical Research/MIT Center for Genome Research, Stanford Human Genome Center, GENETHON 
Human Genome Research Center, and the Cooperative Human Linkage Center. Radiation hybrid (RH) 
mapping was used extensively in the early phase of this study to resolve discrepancies in marker order 
between maps. Specifically, the 83 Stanford G3 radiation hybrid panel was used to map all genetic and STS 
markers available from public database as well as those developed specifically for the project. In addition to 
RH mapping, STS-content mapping using BAC (Bacterial Artificial Chromosome) clones from the region of 
interest was also used routinely to determine the marker order and to complete the BAC contig. 

BAC library screening, end sequencing and contig building 

1 1 . Microsatellite and STS markers obtained from public databases were used to screen the 
human BAC library from Research Genetics (HuntsviUe, AL) by PCR or to the BAC library from Genome 
Systems (St. Louis, MO) screen by hybridization according to manufacturers' protocols. BAC DNA from 
positive clones was prepared, and sequences of the BAC ends were obtained by cycle sequencing the BAC 
DNA directly with vector primers T7 and SP6, respectively, PCR primers were designed from non-repetitive 
end sequences and used as STS markers to improve the physical map and the BAC contig construction. The 
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outlying markers from each side of the contigs were used to screen for overlapping BAC clones to extend the 
contigs. 

Construction of randomly sheared libraries from BACs 

12. BAC DNA was sheared to small fragments of desired size range using a nebulizer. After 
shearing, the libraries were constructed using estabUshed techniques. 

Microsatellite and SNP marker development and genotyping 

1 3 . Microsatellite markers were generated by hybridizing oligonucleotide probes for di, tri, and 
tetranucleotide repeats to randomly sheared sub-libraries made from BAC clones using Quicklite non-isotopic 
enzyme induced chemiluminescent reagents from Lifecodes Corp. (Stamford CT) following the 
manufacturer's instructions. Positive clones were sequenced to identify microsatellite sequences and primers 
were then designed from flanking unique DNA sequence. Primers for ampHfying STS markers were also 
designed using BAC end sequences, and random sequences available within the candidate interval when 
extensive sequencing of the randomly sheared libraries were done. Primer sequences are publicly available at 
PNAS Online. 

14. We genotyped the 4 new microsatellites identified by us in sequencing the region. Primer 
sequences are available on request. Genotyping procedures for the microsatellites were performed using 
established techniques. 

15. Single nucleotide polymorphisms (SNPs) were identified using SSCP (Single Strand 
Conformational Polymorphism) analysis of STS markers (all < 300 basepairs in length), using established 
techniques. We used four unrelated individuals to screen for each SNP. We genotyped the SNPs in patient 
and control samples using standard SSCP procedures. 

Sequencing of the candidate interval and identification of the candidate genes 

16. In the interval of < 3 cM, located within the SAVA5 -ga203 interval, randomly sheared 
libraries prepared from BACs covering this region were sequenced at lOX coverage to discover all sequence 
information and identify all genes within the interval. More than 10,000 individual sequences from the region 
were compared by BLAST20 with sequences from publicly available databases and were analyzed using 
GRAIL21 to identify potential coding sequences. In addition, sequences were assembled using PHRAP 22, 
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23, 24 in a single DNA strand of -33 1 Kb. The whole sequence was again analyzed using BLAST and 
GRAIL to aid in gene prediction. These data were displayed in ACEdb (data available from 
ncbi.nlm.nih.gov) to visualize predicted exons and their relationships to each other. 

Statistical analyses 

17. We applied a modified version of Terwilliger's likelihood ratio test of LD to the 4 novel 
microsateUites and 26 SNPs that spanned our 33 1 Kb candidate region. For each of these 30 markers we applied 
this test twice, once in the sample of 227 patients and their available relatives, and also with the addition of the 
independent controls to the 227 patients and relatives. This likelihood ratio test estimates a single parameter, X, 
which quantifies the overrepresentation of an associated marker allele on disease chromosomes versus control 
chromosomes. X is related to the common epidemiological parameter of population attributable risk. If the 
frequency of an associated allele on disease and normal chromosomes is given by po and Pn, respectively, then X 
is calculated by (po - Pn)K^ - Pu)- Only positive associations with disease are permitted, and X ranges from 0 
(xmder the null of no association) to 1.0 (all disease chromosomes carry the associated allele). Others have shown 
that X is the most closely related to the recombination fraction with disease and less influenced by marker allele 
frequencies than other measures of LD. Because we do not know which chromosome of an affected individual 
harbors the disease locus, we incorporated a genetic model of disease transmission in liie procedure of TerwiUiger. 
Using this model also enabled us to employ data from additional family members other than parents, if th^ were 
not available. The same genetic model (mostly dominant with reduced penetrance) was used as in our previous 
LD papers and in the genome screen of the Costa Rican pedigrees described in Mclnnes et al. In this model one 
chromosome of the affected individual is used as a confrol chromosome. The use of a model is likely to increase 
the power of the test and the precision of the estimates of A, when the inheritance pattern is approximately known. 
Using simulated data, TerwiUiger shows that his test is conservative. 

Results 

Marker development and physical map 

1 8 . Based on our previous results (as described in the instant patent application; and in the 
publications provided herewith) we focused marker development and physical mapping efforts (including 
direct sequencing) in the <3 cM region between sAVA5 and D18S 123 1. Within this region we identified 4 
new microsatellite markers and 26 SNPs to add to the 4 pubhcly available microsatellite markers already used 
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(see Exhibit 3). Based on the extent of haplotype sharing in pedigree CROOl and LD results from the 
previously used markers, we focused our detailed investigation on the region of about 33 1 Kb between PH33 
and D 1 8S 123 1 (although in public databases this segment is estimated as being 378 Kb in length, contig 
NT Ol 1005). Using several sequence analysis tools and database mining procedures (see Methods, above), 
we determined that this interval contained six known genes {CENTRIN, CLULl, TYMS, rTS, YESl, and 
ADCYAPl, ordered from telomeric to centromeric, with TYMS and rTS overlapping each other). This order 
differs in the public database (CENTRIN, CLULl, YESl, rTS, TYMS, mA ADCYAPl, with no overlap 
between rTS and TYMS). All of the genes except "clusterin-like 1 (retinal)" gene {CLULl^ have been well 
characterized previously. CLULl was originally identified during a screen of a human retinal cDNA library 
for retina-specific genes. The function of this gene is not known; however Northern blot analysis reveals that 
it is highly expressed in retina with much lower yet detectable expression in several other tissues including 
brain, kidney and testes. 

Genotyping results 

1 9. We genotyped the 30 new markers in pedigree CROO 1 and in the CVCR patient and control 
samples. Results of the LD analysis for these markers (and the four previously available markers reported in 
ref 8) are displayed in Table 1 (provided herewith as Exhibit 4). Of the 34 markers presented in Table 1,16 
showed association (X >0) with BP-I in at least one of the two samples (that with 227 patients/relatives and 
that with 227 patients/relatives and the addition of 52 controls). The p-value associated with the estimate of 
X was <0.01 for five of these 16 markers, and for four of the five markers the magnitude of association was 
greater in the sample containing the population controls. All five of these markers (PH84, PH205, PH202, 
PH208, and TS30), had estimates of X near 1.0, indicating that virtually all affected individuals had at least 
one copy of the associated allele. The markers showing LD are clustered in the 19 Kb segment between exon 
8 of CLULl and exon 1 of TYMS. This segment also contains the minimal region of haplotype sharing 
within CROO 1, and for each marker in this segment, the associated alleles seen in the population samples are 
the same alleles in the shared haplotype in CROOl (last column in Table 1). 
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Summary 



20. The data presented herein extend the findings described in the instant patent appHcation. The 
patent application provided evidence, from both pedigree analyses and population studies, that a number of 
polymorphisms, including a 154 bp allele of the microsatellite marker D18S59 and a 271 bp allele of the 
microsatellite marker D18S476, are associated with BP. The patent application described how to identify 
additional markers, and how to determine whether such markers are associated with BP. 

2 1 . The data presented herein show that, using techniques described in the patent application, several 
new poljTnorphisms, located in the previously identified interval and associated with BP, were identified. 
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Conclusion 

22 . Those in the field, given the guidance in the instant patent application, could identify 
additional polymorphisms associated with BP. 

23 . I hereby declare that all statements made herein of my own knowledge are true and 
that all statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title XVIII of the United States 
Code, and that such will false statements may jeopardize the validity of the application or any patent 
issuing thereon. 



Date 





Enclosures; Exhibits 2-4 
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Assessing the Feasibility of Linkage Disequilibrium Methods for Mapping 
Complex Traits: An Initial Screen for Bipolar Disorder Loci on 
Chromosome 18 

Michael A. Escamilla/'^'*'+L. Alison Mclnnes/'2'*Mitzi Spesny,*^ Victor I. Reus/ 

Susan K. Service/ Norito Shimayoshi/'*David J. Tyier/ Sandra Silva/ Julio Molina/ 

Alvaro Gallegos/'^Luis Meza/ Maria L Cruz/ Steven Batki/ Sophia Vinogradov/ 

Thomas Neylan/ Jasmine B. Nguyen/ Eduardo Fournier/ Carmen Araya/ 

Samuel H. Barondes/ Pedro Leon/ Lodewijk A. Sandkuijl/ *'' and Nelson B. Freimer'-^ 

'Neurogenetics Laboratory, ^Center for Neurobiology and Psychiatry, ^Department of Psychiatry, San Francisco General Hospital, and 
"Department of Psychiatry, Veterans Administration Medical Center, University of California San Francisco, San Francisco; 'Department of 
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Genetics, Croningen University, Groningen, the Netherlands; and "Cell and Molecular Biology Research Center and Escuela de Medicina, 
Universidad de Costa Rica, and ""Hospital Calderon Guardia, San Jose, Costa Rica 



Summary 

Linkage disequilibrium (LD) analysis has been promoted 
as a method of mapping disease genes, particularly in 
isolated populations, but has not yet been used for ge- 
nome-screening studies of complex disorders. We present 
results of a study to investigate the feasibility of LD 
methods for genome screening using a sample of indi- 
viduals affected with severe bipolar mood disorder (BP- 
I), from an isolated population of the Costa Rican cen- 
tral valley. Forty-eight patients with BP-I were genotyped 
for markers spaced at ~6-cM intervals across chromo- 
some 18. Chromosome 18 was chosen because a pre- 
vious genome-screening linkage study of two Costa Ri- 
can families had suggested a BP-I locus on this 
chromosome. Results of the current study suggest that 
LD methods will be useful for mapping BP-I in a larger 
sample. The results also support previously reported 
possible localizations (obtained from a separate collec- 
tion of patients) of BP-I-susceptibility genes at two dis- 
tinct sites on this chromosome. Current limitations of 
LD screening for identifying loci for complex traits are 
discussed, and recommendations are made for future 
research with these methods. 
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Introduction 

Identifying genes for disorders with complex inheritance 
patterns is one of the greatest challenges in biomedical 
research (Lander and Schork 1994). Such disorders, 
which include many of the most prevalent human dis- 
eases, are difficult to map with standard linkage meth- 
ods. If has been suggested that the availability of dense 
marker maps covering the genome will make linkage 
disequilibrium (LD) analysis a feasible approach for 
screening the genome to map complex disorders (Risch 
and Merikangas 1996). Current marker maps are not 
sufficiently dense to enable such studies to be performed 
in heterogeneous populations or in populations that 
were founded in the distant past. However, the success 
of genome-screening LD-mapping studies of genetically 
simple and/or rare diseases in recently founded isolated 
populations (Houwen et al. 1994; Puffenberger et al. 
1994; Friedman et al. 1995; Newport et al. 1996) pro- 
vide the impetus for testing the utility of LD methods 
for mapping complex diseases in. such populations (Es- 
camilla et al. 1996). In populations where randomly 
sampled patients are on average <20 generations re- 
moved from their last common ancestor, LD may be 
maintained for sizable regions around disease genes. 
Such LD should be manifested by affected individuals 
sharing alleles, identical by descent (IBD), at markers 
spaced at intervals of several centimorgans surrounding 

♦ These authors contributed equally to this work. 

'^Present affiliation: Department of Psychiatry, University of Texas 
Health Science Center at San Antonio, San Antonio, Texas. 

♦Present affiliation: Department ' of Psychiatry, Nara Medical 
School, Nara, Japan. 

^Dr. Gallegos died after this paper was accepted for publication. 
This paper is dedicated to his memory. 



1670 



f Escamilla et al.: LD Mapping of BP-I in Costa Rica 

a disease gene. We now present the results from the first 
stage of a study in which LD methods were used to 
screen for loci that predispose to severe bipolar mood 
disorder,(BP-I), which is common and is almost certainly 
characterized by a complex mode of inheritance. The 
study was done in a relatively recently founded isolated 
population, that of the central valley of Costa Rica 
(CVCR) (Escamilla et al. 1996), where founder effects 
have already been observed for several inherited diseases 
(Saborio 1992; Uhrhammer et al. 1995; Shah et al. 
1997). 

Despite long-standing evidence that BP-I has a genetic 
basis (Escamilla et al. 1997), genome scans for linkage 
have provided equivocal results (Risch and Botstein 
1996; Nurnberger et al. 1997) that fail to satisfy the 
levels of significance suggested for genomewide screens 
by Lander and Kruglyak (1995). The failure to identify 
BP-I loci definitively, by standard linkage approaches, 
probably reflects uncertainty regarding mode of inheri- 
tance, high phenocopy rates, difficulty in demarcation 
of distinct phenotypes, and presumed genetic heteroge- 
neit}'. LD-based mapping approaches within population 
isolates may offer a means of diminishing several of these 
obstacles. An approach (such as LD mapping) that sam- 
ples individuals from an entire population can more eas- 
ily ascertain a large set of patients with a narrowly de- 
fined, reliably diagnosed phenotype (in this case, BP-I) 
than hnkage-based approaches that require ascertain- 
ment of family units with multiple affected cases. Within 
a population isolate, genetic heterogeneity of BP-I may 
also be less than in larger, genetically mixed populations, 
as there is a high probability that individuals with such 
a phenotype share descent from a few common 
ancestors. 

We collected a sample of patients with BP-I, for LD 
analysis, by identifying individuals currently living in the 
CVCR who had known CVCR ancestry. This sample 
was collected independently of our previous pedigree- 
based studies of BP-I in Costa Rica. Our aim in the 
current study was to evaluate the feasibility of identi- 
fying BP-I loci by LD screening in this population, as 
proposed in Escamilla et al. (1996). To do this, we con- 
ducted an LD screen of an entire single chromosome 
(chromosome 18). This chromosome was chosen be- 
cause previous linkage studies in Costa Rica and in other 
populations suggested that it possibly contained bipolar 
disorder loci (Berretini et al. 1994; Stine et al. 1995; 
Freimer et al. 1996a). Genealogical studies indicated that 
the individuals in our current study did not share com- 
mon ancestry over the past several generations (Escam- 
illa et al. 1996). We therefore anticipated that we would 
not detect random genome regions shared IBD by more 
than a few individuals and that regions of high IBD 
sharing would thus be areas containing possible BP- 
I-susceptibility genes inherited from a common founder. 



Samples and Methods 

Sample Collection 

To diminish the likelihood of investigating phenocop- 
ies, we limited the sample to individuals with a definite 
diagnosis of BP-I, with onset by age 50 years and a 
history of at least two psychiatric hospitalizations. The 
48 patients with BP-I (25 female patients and 23 male 
patients) in the current study were recruited indepen- 
dently from psychiatric hospitals and clinics in the 
CVCR. First-degree relatives of patients were also re- 
cruited, to determine genetic phase. The study was ap- 
proved by institutional review boards at the Costa Rican 
Ministry of Health, the University of Costa Rica, and 
the University of California at San Francisco, and in- 
formed consent was obtained from all participating sub- 
jects. Of the 48 BP-I subjects, 8 individuals had both 
parents available for genotyping, 20 individuals had one 
parent available, 10 individuals had one or more chil- 
dren available, 1 individual had two siblings available, 
and 9 individuals had no relatives available. In nuclear 
families, only one individual (the proband) was desig- 
nated as affected, and all others were considered to have 
unknown phenotype. Details of ascertainment and di- 
agnostic procedures, and the cfinical and genealogic pro- 
files of the study sample, can be found in Escamilla et 
al. (1996). 

Genotyping 

We used 26 markers, spanning chromosome 18, to 
genotype all 48 affected individuals (as well as 53 rel- 
atives, to establish phase). Of the 25 regions, 21 were 
^6 cM, and 4 were 6-7 cM. The average distance be- 
tween markers was 4.8 cM. When choices were avail- 
able, we chose the most polymorphic marker (Gyapay 
et al. 1994). The average heterozygosity of the markers 
used in this screen (in the CEPH pedigree collection) was 
0.75. (The only screening markers with heterozygosity 
values <0.70 were D18S464, D18S60, D18S378, and 
D18S469.) We screened chromosome 18 at a marker 
density of 6 cM because available marker maps had gaps 
^6 cM, and our goal was to have an equal density of 
coverage across the chromosome (Gyapay et al. 1994). 
We chose markers from the maps available, at the time 
of the study, from Genethon (Gyapay et al. 1994), the 
Cooperative Human Linkage Center (Murray et al. 
1994), and the public database of the Utah Center for 
Genome Research. Genotyping procedures used for all 
experiments were as previously described by Di Rienzo 
et al. (1994). In brief, one of the two primers was labeled 
radioactively with a polynucleotide kinase, and PCR 
products were separated, by electrophoresis, onto poly- 
acrylamide gels. Autoradiographs were scored indepen- 
dently by two raters. Data for each marker were entered 
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into the computer database twice, and the resultant files 
were compared for discrepancies. Scoring was done 
without knowledge of affected status. 

Simulations 

We conducted simulations to evaluate the power of a 
likelihood-based test of LD (Terwilliger 1995), to detect 
a result significant at the .05 level, with these assump- 
tions: a 6-cM marker map; a disease gene in the middle 
of the 6-cM segment; affected subjects, with one copy 
of the disease gene, separated by 10 generations from a 
common ancestor; and four equally frequent marker al- 
leles at each marker site. (The disease gene was associ- 
ated with the "1 " allele at the marker locus.) Under these 
assumptions, and with a phenocopy rate of 0%, normal 
chromosomes carried each marker allele with a proba- 
bility of 25% (normal-chromosome distribution), and 
disease chromosomes carried the " 1 " allele with a prob- 
ability of 80%. The probability of disease-chromosome 
distribution was calculated with the formula (1 — 
ef -t- [1 - (1 - df X f], where 6 = recombination 
fraction, G = number of generations from a common 
ancestor, and f= the frequency of the allele in the pop- 
ulation. Thus, the disease chromosomes carry the "1" 
allele with a probability of 80% and each of the re- 
maining three alleles with a probability of 6.7%. Because 
the true genetic structure of bipolar disorder is unknown, 
we examined several different conditions of etiologic het- 
erogeneity (which would include locus and allelic het- 
erogeneity, as well as phenocopies). We investigated phe- 
nocopy rates of 0%, 33%, and 67% (with phenocopy 
rates of 33% and 67%, the percentages of chromosomes 
from affecteds with the "1" allele are 62% and 43%, 
respectively). If an affected individual was randomly se- 
lected as a phenocopy (with a probability equal to the 
phenocopy rate), then the marker allele on all four pa- 
rental chromosomes was randomly chosen from the nor- 
mal chromosome distribution. If the affected individual 
was randomly chosen as a true case, (with a probability 
of 1 minus the phenocopy rate) the marker allele for 
one chromosome of that individual was randomly cho- 
sen from the normal chromosome distribution, and the 
other chromosome's marker alleles were randomly cho- 
sen from the disease-chromosome distribution. Recom- 
bination occurred on parental chromosomes in propor- 
tion to the marker map. Marker alleles for 
nontransmitted chromosomes of the parents were ran- 
domly chosen from the normal chromosome distribu- 
tion. We performed these analyses by using the 48 pa- 
tients with BP-I plus their available relatives. One 
hundred replications were performed for each simula- 
tion. Available relatives were considered to have un- 
known disease phenotype. For the 10 affected individ- 
uals with at least one child available for genotyping, one 
chromosome from the affected parent was randomly 



simulated to be transmitted to available children, and 
the other chromosome was randomly selected from the 
normal chromosome distribution. Although data were 
simulated for parents of all affected individuals, if par- 
ents were not available for genotyping, their simulated 
genotypes were not used in these analyses. 

We also did power simulations (100 replications for 
each model) of larger sample sizes, using an ideal situ- 
ation in which both parents are available for genotyping, 
to aid in planning future studies. In these simulations 
we used sample sizes of 90, 200, 300, and 400 affected 
individuals; phenocopy rates of 50% and 75%; and a 
marker map of 2.5 cM, with all other assumptions as 
described above. With this denser marker map, at a phe- 
nocopy rate of 0%, disease chromosomes carried the 
"1" allele with a probability of 90%, calculated by the 
formula {1 - Of + [1 - {1 - ef x f], and each of the 
remaining three alleles with a probability of 3.3%. De- 
tails of the likelihood-ratio test used in analyzing sim- 
ulation results are described inAnalysis. 

Analysis 

We used two different procedures to identify regions 
potentially shared IBD by patients with BP-I. The first 
approach, a search for shared segments, has the advan- 
tage of being nonparametric. The second approach, al- 
though requiring parameters of the illness to be specified, 
has the advantage of providing a formal test statistic, 
allowing for the calculation of P values. These two tests 
thus offer compensatory strengths and weaknesses when 
used in the search for genes in a complex disease. 

We first searched for shared segments (Houwen et al. 
1994). For each individual, we evaluated two marker 
haplotypes in each of the 25 intermarker intervals, by 
using a preselected threshold (the possible sharing of a 
haplotype by ^50% of patients) to select segments for 
■further investigation. Since this screen does not differ- 
entiate between sharing that is IBD and sharing that is 
identical by state (IBS), use of lower thresholds would 
lead to too many segments passing the screen. 

We also applied a likelihood-ratio test for LD to each 
of the 26 initially tested markers. This test was done 
independently of the results of the shared-segment eval- 
uation. We applied a modified version of the procedure 
of Terwilliger (1995), which only includes case and con- 
trol chromosomes or chromosomes transmitted and not 
transmitted to patients. In our sample there were several 
affected individuals whose parents were not available 
but whose children were available. DNA from these lat- 
ter individuals could not be analyzed with the original 
Terwilliger program but could be analyzed with our im- 
plementation of the same procedure, as described by 
Freimer et ai. (1996a). This procedure examines the like- 
lihood that a particular allele (or alleles) is (are) over- 
represented on disease chromosomes compared with 
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Table 1 

Hoterozygosity of Markers Used in the Genome Screen of 
Chromosoine 18 



Marker Name 


Genethon Database 


Costa Rican Sample* 


_ 

1)]SS1140 


.49 


.39 


D18S59'' 


.81 


.81 


D18S476'' 


.76 


.62 


D]SS481 


...76 


.74 


D18S391 


.75 


.69 


D18S452 


.83 


.85 


D18S843 


NA 


.73 


DISS464 


.65 


.51 


D18S1153 


.78 


.69 


1)18S378 


NA 


.54 


D18S53 


.79 


.81 


D18S453 


.82 . 


.81 


D18S40 


NA 


.81 


D18S66 


.85 


.81 


D18S56 


.73 


.74 


D18S57 


.87 


.85 


0185467" 


.73 


.64 


D18S460 


.62 


.67 


D18S450 


.79 


.74 


D18S474 


.82 


.73 


D18S69 


.79 


.78 


D18S64 


.74 


.65 


D18S1134 


.73 


.68 


D18S1147 


.85 


.86 


D18S60 


.37 


.58 


D18S55 


.77 


.80 


D18S68 






D18S477 


.62 


.70 


Dissei" 


.87 


.86 


D18S488 


.87 


.82 


DI8S485'' 


.79 


.79 


D18S541 


NA 


.63 


D18S870'' 


NA 


.66 


D18S469'' 


.65 


.64 


D18S874 


NA 


.64 


D18S380 


NA 


.63 


01851121" 


.74 


.77 


D18S1009 


.74 


.66 


D18S844 


NA 


.76 


DI8SJ54 


.82 


.79 


D18S461 


.77 


.65 


D18S70 


.83 


.86 



Note. — NA = data not available. ■ 

' Allele frequencies were calculated from the entire sample, a 
:ounting for known relationships among individuals. 
Markers with -21«(LR) >1.0. 
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is nearly dominant (assuming penetrance of .81 for het- 
erozygous individuals and .90 for homozygous individ- 
uals with the disease mutation), that the population 
prevalence of BP-I in Costa Rica is .015, and that the 
frequency of the disease gene in the population is .003. 
In the likelihood calculations, all possible disease-ge- 
notype combinations of all relatives are considered. "With 
the model that was used, in which the disease-gene fre- 
quency is very low, the LD-likelihood test, in most cases, 
treats the probands as effectively heterozygous at the 
disease locus, and chromosomes of other relatives not 
occurring in the probands are treated as non-disease- 
carrying chromosomes. We did not specify a phenocopy 
rate in the genetic model, because the effect of pheno- 
copies will be absorbed by the parameter X; the presence 
of phenocopies in our sample will serve to erode the 
association between marker alleles and disease and 
hence will reduce the estimate of X. Because, in the pre- 
sent LD study, we were attempting to gather further 
evidence regarding the findings published in our initial 
genome screen, we limited ourselves to this one model 
in performing the likelihood analyses. However, both 
the BP-I family sample and the current LD sample will 
ultimately be analyzed with use of other models. "We 
considered as promising those markers that gave evi- 
dence of overrepresentation of an allele on affected chro- 
mosomes, with a — 21w(likelihood ratio [LR]) statistic 
>1.0. 

Follow-up genotyping and LD-analysis studies were 
performed on markers that gave suggestive findings in 
the shared-segment evaluation. "Within each segment 
that passed the threshold described above, 1-3 addi- 
tional markers were typed to permit us to test for LD 
across regions of 1-2 cM. Markers that provided sug- 
gestive evidence of LD by the initial likelihood-ratio test, 
but had- not been suggested as promising by the shared- 
segment Screen, were also followed up, in this case by 
typing two additional nearby markers. In all, a total of 
42 markers from chromosome 1 8 were used to genotype 
the study sample (table 1 and fig. 1). LD analysis of the 
additionally typed markers was conducted by use of the 
likelihood-ratio test. 

Results 



nondisease chromosomes (Terwilliger 1995; Freimer et 
al. 1996<3). A single parameter, X, is estimated, which 
quantifies such overrepresentation of marker alleles on 
I disease chromosomes. Designation of chromosomes of 
f probands as disease carrying or non-disease carrying 
I was achieved by specification of a genetic model for the 
1 disease. The same model of transmission was used in 
C this LD-likelihood test as was used in the initial genome 
I screen of the Costa Rican families, described in Mclnnes 
^ et al. ( 1 996). In brief, this model assumes that the disease 

k 



Simulations 

Simulation results for the sample of 48 patients with 
BP-I and available relatives showed relatively high power 
to detect suggestions of association (P ^ .05) with low 
phenocopy rates (94% for a phenocopy rate of 0%, and 
54% for a phenocopy rate of 33%) but a dramatically 
decreased power under high phenocopy rates (e.g., 9% 
for a phenocopy rate of 67%). Additional simulations 
showed that, under higher phenocopy rates, the power 
to detect LD can be improved by increasing the sample 



1674 



Am.]. Hum. Genet. 64:1670-1678, 1999 



-2*Ln(LikGiihood Ratio) 



pter 




Figure 1 Results from the LD screen of chromosome 18. The 
26 markers used in the first stage of the screen are listed in the right 
column. Sixteen markers used to follow up interesting regions are listed 
in the left column. Approximate chromosomal locations of the 26 
initial markers and the 16 follow-up markers are indicated by long 
and short tick marks, respectively. The eight segments that passed the 
initial screen threshold for segment sharing (50% of individuals or 
25% of chromosomes sharing a two-marker haplotype) and the five 
markers that passed the initial threshold for the Terwilliger LR test 
(-2l«[LR] > 1.0) are indicated by blackened bars and asterisks, re- 
spectively. Two marker segments that passed the initial threshold were 
followed up by at least one marker within the segment, if possible (at 
the time of the study no markers were available between D18S843 
and D18S464, and only one marker was available between D18S464 
and D18S378). Markers that passed the initial threshold for the Ter- 
williger LR test were followed up with two additional markers. These 
additional markers flanked the original finding. The value of the 
-2hi(LR) statistic, from the Terwilliger test, is plotted as a solid bar. 
This statistic is distributed as a one-sided random variable with one 
degree of freedom. The estimate of the X value, for the eight markers 
with positive resuhs, is indicated in parentheses after the -21n(LR) 
statistic. Markers without a -21«{LR), statistic plotted had estimates 
of X = 0, with the exception of three markers that had estimates of 
0<X<0.62. 

size and/or the marker density of screening (table 2). For 
instance, with a phenocopy rate of 75%, the power in- 
creases to 82% with a sample of 300 affected individuals 
and a 2.5-cM marker map. 



Shared-Segment Screen 

We evaluated 25 possible shared segments (defined by 
the 26 markers genotyped in the sample). Eight regions 
passed the threshold of possible IBD sharing by 3=50% 
of patients. These regions were bounded by the following 
markers: D18S843-D18S464, D18S464-D18S378, 
D18S467-D18S474, D18S64-D18S60, D18S60- 
D18S68,D18S485-D18S469,D18S469-D18S1009,and 
D18S1009-D18S461 (fig. 1). 

Linkage-Disequilibrium Testing 

Five (D18S59, D18S467, D18S61, D18S485, and 
D18S469) of the original 26 markers displayed evidence 
of possible LD, by means of a likelihood procedure 
(-21m[LR] statistic >1.0; table 3). Two (D18S59 and 
D18S61) of these five markers had not been identified 
as markers of interest by the shared-segment evaluation. 
D18S59, located near ISpter, displayed the strongest 
pointwise evidence for LD (-2ln[LR] statistic of 8.3, 
P = .002) of all the markers tested in this sample. 

Follow-up of Initial Results 

Using the protocol discussed in Samples and Methods, 
we genotyped additional markers within the segments 
that passed the shared-segment screen as well as follow- 
up markers surrounding one (D18S59) of the two mark- 
ers that had passed only the LD screen. We were unable 
to follow up one shared-segment region (D18S843- 
D18S464), because additional polymorphic markers 
were not available within the segment. We were also 
unable to follow up the finding for Dl 8S61, for the same 
reason. Three (D18S476, D18S870, and D18S1121) of 
the 16 follow-up markers typed displayed additional ev- 
idence of possible LD (fig. 1). 

These additional results brought to eight the total 
number of markers with -21n(LR) statistics >1.0 (table 
3). Five of these eight marker loci were clustered within 
a small region of 18q22-23. The most significant LD in 
18q22-23 was observed at D18S1121, with -2ln(LR) 
of 5.03 and P = .01, and two were in ISpter. 

For the two ISpter markers (D18S59 and D18S476), 

Table 2 

Power-of-Likelihood-Analysis Test of LD 

Power to detect LD for 
Sample Size (n) = 

Penocopy Rate 90 ■ 200 300 400 

50% ^ -81% 99% 100% 100% 

75% 33% 62% 82% 90% 

Note.— Assumptions included that subjects were removed from a 
common ancestor by 10 generations, that a marker map of 2.5 cM 
was used, and that each marker had four equally frequent alleles. 
Values are the percentage of replicates to have P values <.05. 
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Table 3 

Frequencies of Marker Alleles Overrepresented in Disease 
Chromosomes, as Compared with Nondisease Chromosomes, for 
Markers Where -2ln(LR) >1 .0 

Frequency on 



Nondisease Disease 
Marker Allele Chromosomes Chromosomes 



D18S59' 154 .121 .572 

D18S476 271 - .470 .771 

D18S467" 172 .384 .693 

DlSSei" 177 .074 .326 

D18S485' 182 .237 .586 

D18S870 179 .405 .657 

D18S469' 234 .128 .450 

D18S1121 168 .171 .553 



" Markers from the screening stage. 

the alleles overrepresented on BP-I chromosomes (154 
and 271 bp, respectively) form a haplotype that occurs 
in 48% of the patients with BP-I. Overall, this haplotype 
occurs on 26% of the chromosomes of individuals with 
BP-I and on 4% of the chromosomes not transmitted 
from parents to individuals with BP-I (definite phase for 
these two markers could be assigned in 25 patients with 
BP-I [50 chromosomes] and 25 nontransmitted parental 
chromosomes). Because the composite genetic and phys- 
ical maps of the 18q22-23 region had not yet been com- 
pleted at the time of this study, the relative order of the 
five markers in 18q22-23, for which evidence of LD was 
observed, was still too uncertain to permit construction 
of definitive marker haplotypes in our study sample. 

Marker D18S467, in the 18ql23 region, was the one 
marker outside 18q22-23 and ISpter to show a 
-21n(LR) >1 (-21n[LR] = 2.5, P = .06). The additional 
markers used to follow up this result (D18S450, 
D18S460, and D18S57) displayed no evidence of 
association. 

Marker Heterozygosity in the Costa Rican Sample 

We calculated heterozygosity values for the markers 
used, on the basis of the allele frequencies, estimated 
from the entire sample, accounting for known relation- 
ships among individuals. These heterozygosities are 
shown in table 1, along with the corresponding hetero- 
zygosity values of these markers in the CEPH popula- 
tion, used by Genethon. 

Discussion 

Screening for Complex Disease Loci by LP 
Approaches 

Our intention in this work was to explore the feasi- 
bility of using LD methods to screen the genome for 
susceptibility genes for a common, genetically complex 



disorder. The results obtained in our LD-based search 
for possible BP-I gene-loci on chromosome 18 were en- 
couraging (specific susceptibility regions were sug- 
gested), but they highlight a number of issues that must 
be considered before LD screening is widely adopted. 

Successful application of a shared-segment approach 
to any LD study depends on ( 1 ) a marker-map density 
that is appropriate to the age of the population isolate 
being studied and (2) a sharing threshold that will not 
be too high to allow true IBD areas to be identified and 
that will not be so low as to include many areas that 
are IBS false-positive signals. An appropriate marker 
map for an LD-screening study should have segments of 
a size expected to be shared IBD by many of the affected 
individuals. In addition to the density of the marker map 
used, the number of generations separating affected in- 
dividuals from their common ancestor and the rate of 
etiologic heterogeneity in the population will also influ- 
ence the choice of the sharing threshold, used to trigger 
further study. For example, if the common (disease-gene 
bearing) ancestor is removed from the current descen- 
dants by >10 generations, the length of true IBD hap- 
lotypes shared by 2*50% of the descendants may be <5 
cM (and certainly <6 cM, as is the screen used in this 
study) (Te Meerman et al. 1994; Durham and Feingold 
1997). Our choice of a threshold of 50% of affected 
individuals sharing a possible haplotype therefore effec- 
tively meant that we were likely to identify only BP-I 
genes of a major effect in this population (phenocopy 
rate approaching zero), and even then, only if the dis- 
tance from a common ancestor is not > ~10 generations. 
Although this was probably too stringent a screen 
threshold, given the complex etiology of bipolar disor- 
der, the alternative we faced — reducing the threshold to 
a lower percentage of potential IBD sharing — would 
have drastically decreased the specificity (and hence the 
utility) of the screen. For instance, in this particular 
study, lowering the threshold to a possible IBD haplo- 
type shared in 5=25% of the patients would have resulted 
in 24 of the 25 regions tested being determined as regions 
of interest. If, in future studies, definite phase infor- 
mation can be set for a greater proportion of the pro- 
bands (obtained from phasing information supplied by 
additional relatives) the "possible IBD" threshold will 
be more useful as a screening criteriori at thresholds ap- 
proaching 25% sharing (almost one in five of the pa- 
tients with BP-I in the current study had no relative 
available for phase construction). Finally, regardless of 
the threshold chosen, there is no widely accepted statis- 
tical test available to evaluate the significance of the 
number of shared haplotypes observed, although several 
statistical approaches are under development (reviewed 
by Kruglyak 1997). 

The use of markers with low heterozygosity will in- 
crease the number of false-positive results in a shared- 
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The third region that showed possible evidence of LD 
in our population sample was identified through a single 
marker {D18S467), in the 18ql2.3 region. Additional 
markers typed near this one did not support the initial 
suggestion of LD in this region. Evidence from a linkage 
test that yielded a significance level of ? = .06 would be 
expected to occur, by chance, ~24 times (about once on 
most chromosomes) in a genomewide screen. 

Our possible BP-I localizations at ISpter and 18q22- 
23, in the current sample, are distinct from regions on 
chromosome 1 8 suggested by other groups as being pos- 
sibly linked to mood disorder (Berretini et al. 1994; Stine 
et al. 1995). We detected no evidence of association with 
these areas (near the centromere and in 18q21) in our 
BP-I population sample; nevertheless, the power of our 
current sample is not great enough to rule out these 
regions as potential BP-I loci in the Costa Rican pop- 
ulation. McMahon et al. (1997) have recently reported 
excess allele sharing in sib pairs at 18S541, which is in 
the 1 8q22-23 region, although their affected status in- 
cluded not only BP-I, but also bipolar type II and schi- 
zoaffective patients. 

Future Directions 

The results of this study suggest that shared-seg- 
ment-screening approaches will only be useful with the 
development of denser marker maps (ColHns et al. 1997) 
and with the development of tests that permit statistical 
comparison of disease-chromosome haplotypes with 
control-chromosome haplotypes. Because the potential 
advantages of a shared-segment approach are substantial 
(this type of approach takes maximum advantage of the 
fact that haplotypes, not just single alleles, are inherited 
IBD in population isolates, and it is nonparametric), and 
because marker maps (Dib et al. 1996; Yuan et al. 1997) 
and statistical methods continue to improve, we remain 
optimistic about this method of mapping genes for com- 
plex disorders. 

Both the ISpter and the 18q22-23 regions would have 
been identified as regions with possible LD at a signif- 
icance of P< .05, even if we had not used the shared- 
segment approach but had instead screened for evidence 
by using only the likelihood-ratio test, with the original 
26 markers. Our results indicate that, in population iso- 
lates, such as the CRCV, and with suitably dense marker 
coverage, tests similar to the likelihood-ratio test of LD 
(Terwilliger 1995) are promising tools for genome 
screening of complex diseases. It is not clear, however, 
whether currently available tests will be powerful 
enough to detect unequivocal proof of association in a 
genomewide scan for such diseases, given sample sizes 
that are easily obtained. More-powerful tests are needed 
and may emerge from efforts to develop measures that 
make use of haplotype information (Service et al. 1999 



[in this issue]; Durham and Feingold 1997; Goldin and 
Chase 1997). 

The test of LD screening conducted in the current 
study points out the need to do a more complete LD 
screening analysis. We thus intend to perform an LD 
screen of chromosome 18, using an expanded sample of 
patients with BP-I and a denser marker map. The ad- 
dition of more-polymorphic markers to genome maps, 
and the application of haplotype-based statistical tests 
currently under development, should facilitate efforts to 
definitively identify BP-I susceptibility genes in Costa 
Rica. 

Acknowledgments 

This work was supported by National Institutes of Health 
grants MH00916, MH49499, MH48695, and MH47563; the 
National Alliance for Research on Schizophrenia and 
Depression Young Investigator Awards (to M.A.E. and N.B.F.); 
the American Psychiatric Association's Program for Minority 
Research Training in Psychiatry Award (to M.A.E.); and the 
Veterans Administration Research Psychiatrist Award (to 
L.A.M.). We thank Lauren Smith for assistance in collating 
the best estimate materials. M.A.E. thanks Dr. Carlos Martinez 
for his support'and suggestions regarding this research work. 
Above all we would like to thank the patients and family 
members who participated in this project and the Costa Rica 
institutions that made this work possible, including the Hos- 
pital Nacional Psiquiatrico, The Hospital Calderon Guardia, 
The Ministry of Science, and the Caja Costarricense de Seguro 
Social. 

Electronic-Database Information 

URLs for data in this article are as follows: 

Cooperative Human Linkage Center, http://www.chlc.org/ (for 

marker maps) 
Genethon, http://www.genethon.fr/ (for marker maps) 
Utah Center for Genome Research, http://www.genome.utah 

.edu/ (for marker maps) 

References 

Berretini WH, Ferraro TN, Goldin LR, Weeks DE, Detera- 
Wadleigh S, Nurnberger JI Jr, Gershon E (1994) Chromo- 
some 18 DNA markers and manic-depressive illness: evi- 
dence for a susceptibilky gene. Proc Natl Acad Sci USA 91: 
5918-5921 

Collins FS, Guyer MS, Chakravard A (1997) Variations on a 
theme: cataloging human DNA sequence variation. Scierice 
278:1580-1581 . ■ ' 

Dib C, Faure S, Fizames C, Samson D, Prouot N, Vignal A, 
Millasseau P, et al (1996) A comprehensive genetic map of 
the human genome based on 5,264 microsatellites. Nature 
380:152-154 

Di Rienzo A, Peterson AC, Garza JC, Valdes AM, Slatkin M, 
Freimer NB (1994) Mutational processes of simple-sequence 



1678 



Am. J. Hum. Genet 64:1670-1678, 1999 



repeat loci in human populations. Proc Natl Acad Sci USA 
91:3166-3170 

Durham LK, Feingoid E (1997) Genome scanning for segments 
shared identical by descent among distant relatives in iso- 
lated populations. Am J Hum Genet 61:830-842 

Escamilla M, Spesny M, Reus V, Gailegos A, Meza L, Molina 
J, Sandkuijl L, et ai (1996) Use of linkage disequilibrium 
approaches to map genes for bipolar disorder in the Costa 
Rican population. Am J Med Genet 67:244-253 

Escamilla MA, Freimer NB, Reus VI (1997) The genetics of 
bipolar disorder and schizophrenia. In: Rosenberg R, Pru- 
siner S (eds) Molecular and genetic basis of neurological 
disease, 2d ed. Butterworth Heinemann, Newton, pp. 
1343-1362 

Friedman TB, Liang Y, Weber JL, Hinnant JT, Barber TD, 
Winata S, Arhya IN, et al (1995) A gene for congenital, 
recessive deafness DFNB3 maps to the pericentromeric re- 
gion of chromosome 17. Nat Genet 9:86-91 

Freimer NB, Reus VI, Escamilla MA, Mclnnes LA, Spesny M, 
Leon P, Service SK, et al {\996a) Genetic mapping using 
haplotype, association and linkage methods suggests a locus 
for severe bipolar disorder (BPI) at 18q22-q23. Nat Genet 
12:436-441 

Freimer NB, Reus VI, Escamilla M, Spesny M, Service S, Ga- 
ilegos A, Meza L, et al (1996b) An approach to investigating 
linkage for bipolar disorder using large Costa Rican pedi- 
grees. Am J Med Genet 67:254-263 

Goldin LR, Chase GA (1997) Improvement of power to detect 
complex disease genes by regional inference procedures. Ge- 
net Epidemiol 14:785-789 

Gyapay G, Morisette J, Vignal A, Dib C, Fizames C, Millasseau 
P, Marc S, et al (1994) The 1993-1994 Genethon human 
genetic linkage map. Nat Genet 7:246-339 

Houwen RHJ, Baharloo S, Blankenship K, Raeymakers P, Juyn 
J, Sandkuyl LA, Freimer NB (1994) Genome screening by 
searching for shared segments; mapping a gene for benign 
recurrent intrahepatic cholestasis. Nat Genet 8:380-386 

Kruglyak L (1997) What is significant in whole-genome link- 
age disequilibrium studies? Am J Hum Genet 61:810-812 

Lander E, Kruglyak L (1995) Genetic dissection of complex 
traits: guidelines for interpreting and reporting linkage re- 
sults. Nat Genet 11:241-247 

Lander E, Schork NJ (1994) Genetic dissection of complex 
traits. Science 265:2037-2948 

Mclnnes LA, Escamilla MA, Service SK, Reus VI, Leon PE, 
Silva S, Rojas E, et al (1996) A complete genome screen for 
genes predisposing to severe bipolar disorder in two Costa 
Rican pedigrees. Proc Natl Acad Sci USA 93:13060-13065 

McMahon FJ, Hopkins PJ, Xu J, Mclnnes MG, Shaw S, Car- 
don L, Simpson SG, et al (1997) Linkage of bipolar affective 
disorder to chromosome 18 markers in a new pedigree series. 
Am J Hum Genet 61:1397-1404 

Murray JC, Buetow KH, Weber JL, Ludwigsen S, Scherpbier- 
Heddema T, Manion F, QuiUen J, et al (1994) A compre- 



hensive human linkage map with centimorgan density. Sci- 
ence 265:2049-2054 

Newport MJ, Huxley CM, Huston S, Hawrylowicz CM, Oos- 
tra BA, Wilhamson R, Levin M (1996) A mutation in the 
interferon-7-receptor gene and susceptibility to mycobac- 
terial infection. New Engl J Med 335:1941-1949 

Nurnberger JI, DePaulo RJ, Gershon ES, Reich T, Blehar MC, 
Edenberg HJ, Foroud T, et al (1997) Genomic survey of 
bipolar illness in the NIMH genetics initiative pedigrees; a 
preliminary report. Am J Med Genet 74:227-237 

Puffenberger EG, Kauffman ER, Bold S, Matise TC, Wash- 
ington SS, Angrist M, Weissenbach J, et al (1994) Identity- 
by-descent and association mapping of a recessive gene for 
Hirschsprung disease on human chromosome 13q22. Hum 
Mol Genet 3:1217-1225 

Risch N, Botstein D (1996) A manic depressive history. Nat 
Genet 12:351-353 

Risch N, Merikangas K (1996) The future of genetic studies 
of complex human diseases. Science 273:1516-1517 

Saborio M (1992) Experience in providing genetic services in 
Costa Rica. In: Birth Defects Original Article Series 28: 
96-102 

Service SK, Temple Lang DW, Freimer NB, Sandkuijl LA 
(1999) Linkage-disequilibrium mapping of disease genes by 
reconstruction of ancestral haploytpes in founder popula- 
tions. Am J Hum Genet 64:1729-1739 (in this issue) 

Shah AB, Chernov I, Zhang HT, Ross BM, Das K, Lutsenko 
S, Parano E, et al (1997) Identification and analysis of mu- 
tations in the Wilson disease gene {ATP7B): population fre- 
quencies, genotype-phenotype correlation, and functional 
analysis. Am J Hum Genet 61:317-328 

Spielman RS, McGinnis RE, Ewens WJ (1993) Transmission 
test for linkage disequilibrium: the insulin gene region and 
insulin-dependent diabetes mellitus (IDDM). Am J Hum Ge- 
net 52:506-516 

Stine OC, Xu J, Koskela R, McMahon FJ, Gschwend M, Frid- 
dle C, Clark CD (1995) Evidence for linkage of bipolar 
disorder to chromosome 18 with a parent-of-origin effect. 
Am J Hum Genet 57:1384-1394 

Te Meermah GJ, van der Meulen MA, Sandkuijl LA (1994) 
Expected size of shared haplotypes surrounding a common 
disease gene. Am J Hum Genet Suppl 55;A205 

Terwilliger JD (1995) A powerful hkelihood method for the 
analysis of linkage disequilibrium between trait loci and one 
or more polymorphic marker loci. Am J Hum Genet 56: 
777-787 

Uhrhammer N, Lange E, Porras O, Naeim A, Chen X, Sheik- 
havandi S, Chiplunkar S (1995) Sublocalization of an ataxia- 
telangiectasia gene distal to D11S384' by ancestral haplo- 
typing in Costa Rican families. Am J Hum Genet 57: 
103-111 

Yuan B, Vaske D, Weber JL, Beck J, Sheffield VC (1997) Im- 
proved ' set of short-tandem-repeat polymorphisms for 
screening the human genome. Am J Hum Genet 60:459-460 



